Antibodies utilizing VL6-57 light chains target a convergent cryptic epitope on SARS-CoV-2 spike protein and potentially drive the genesis of Omicron variants

Continued evolution of SARS-CoV-2 generates variants to challenge antibody immunity established by infection and vaccination. A connection between population immunity and genesis of virus variants has long been suggested but its molecular basis remains poorly understood. Here, we identify a class of SARS-CoV-2 neutralizing public antibodies defined by their shared usage of VL6-57 light chains. Although heavy chains of diverse genotypes are utilized, convergent HCDR3 rearrangements have been observed among these public antibodies to cooperate with germline VL6-57 LCDRs to target a convergent epitope defined by RBD residues S371-S373-S375. Antibody repertoire analysis identifies that this class of VL6-57 antibodies is present in SARS-CoV-2-naive individuals and is clonally expanded in most COVID-19 patients. We confirm that Omicron-specific substitutions at S371, S373 and S375 mediate escape of antibodies of the VL6-57 class. These findings support that this class of public antibodies constitutes a potential immune pressure promoting the introduction of S371L/F-S373P-S375F in Omicron variants. The results provide further molecular evidence to support that antigenic evolution of SARS-CoV-2 is driven by antibody mediated population immunity.

The SARS-CoV-2 pandemic offered an opportunity to interrogate immune response towards an emergent virus infecting an immunologically naive population.During the first 3 years of the pandemic, the spike (S) protein on the virus surface has been actively undergoing antigenic drift generating a large number of variants [1][2][3][4] .In late 2021, the Omicron BA.1 variant emerged, its highly mutated S-protein renders most antibodies isolated earlier in the pandemic ineffective, conferring the Omicron BA.1 variant with an extraordinary immune evasion capability 5,6 .Further antigenic drift of S-protein had been observed and numerous Omicron subvariants emerged, including BA.2-BA.5 [7][8][9] .Since late 2022, XBB and BQ.1 subvariant series have emerged based on BA.2 and BA.5 respectively with further antigenic changes on the S-protein 10,11 .As the pandemic progresses, S-protein is drifting antigenically further generating new variants.
The emergence of Omicron variants marked great shifts in S-protein antigenicity and virological behaviors, sparking several theories on origins of Omicron variants 12 .Several recurrent mutations observed in earlier variants became fixed in the Omicron S-protein, most of these are located within the receptor binding domain (RBD) of the S-protein, including K417N, L452R, T478K, E484A and N501Y.Earlier studies have identified that some of these mutations are located in the epitopes recognized by several classes of germline antibodies that are widely induced within the population.K417N is located in the epitope of VH3-53/3-66 encoded class 1 antibodies inhibiting ACE2 binding [13][14][15] ; L452R is located in the epitope of VH1-69 encoded population antibodies [16][17][18] ; E484A is located in the epitope of VH1-2 encoded class 2 population antibodies [19][20][21] .These mutations often confer viruses with considerable antibody evading ability.While L452R and N501Y have also been proposed to enhance or compensate fusogenic or receptor binding activity of S-protein 4,[22][23][24] , association of recurrent mutations to population antibody responses strongly implies that antigenic drift of SARS-CoV-2 S-protein may be driven by immune pressure at population level.Population immunity (or herd immunity) has long been proposed to drive antigenic drift of influenza virus haemagglutinin 25 , despite molecular basis behind such driving force remains poorly characterized.Apart from the above identified recurrent mutations, several RBD mutations were first emerged with the highly mutated Omicron BA.1 S-protein, including S371L, S373P, S375F, N440K, G446S, S477N, Q493R, G496S and Q498R, most of them have been fixed in the S-protein of current circulating viruses.Among these, S371L (S371F in BA.2-BA.5),S373P and S375F are particularly interesting that they are close together in sequence while located relatively far away from the receptor binding site 26 .The S371L/F-S373P-S375F mutations are considered as an important feature of Omicron variants and have been recently shown as an important determinant for the Omicron nasal cell tropism 27 .They have also been shown to modulate S-protein functions 28,29 and antigen presentation 30 .It has been hypothesized that the altered Omicron tissue tropism can also affect disease severity [31][32][33][34] .However, the driving force behind the genesis of the S371L/F-S373P-S375F mutations remains enigmatic.In this study, we report that diverse antibodies utilizing IGLV6-57 (VL6-57) light chains are widely induced within the population and they target a common epitope defined by S371, S373 and S375 in the ancestral SARS-CoV-2 S-protein.We provide evidence to suggest that the identified class of antibodies constitute a potential immune pressure driving the genesis of S371L/F-S373P-S375F mutations in Omicron variants.

Binding of R1-26 antibody to SARS-CoV-2 S-trimer
To further understand neutralization activity of R1-26, we determined cryo-electron microscopy (cryo-EM) structures of S-trimer:R1-26 Fab complexes in two different stoichiometries (3:2, 3:3 S-protomer:R1-26 Fab) (Figs. 2a, S1a) with upper resolutions in the range of 3.19-3.44Å (Fig. S2a-c).An unusual structure showing a head-to-head aggregate of two S-trimers each bound by 3 R1-26 Fabs was also determined at a lower resolution of 5.3 Å (Figs. 2a, S1a, S2a-c), the aggregation is mediated by interactions between Fabs bound to different S-trimers and a similar aggregate was observed for the S-trimer:IgG complex of 6M6 (using IGHV3-9, IGLV3-21 genes) 35 .It is not known whether the observed interactions could lead to aggregation of S-trimers on virus surface.
The detailed R1-26 epitope structure confirms that the epitope is fully buried in the previously determined RBD "down" closed and locked S-trimer structures [36][37][38] , therefore, the epitope is fully cryptic in a S-trimer adopting 3-RBD "down" closed or locked conformations (Fig. S3a-h).In the 1 RBD "up" S-trimer, this epitope on the "up" RBD is also largely obstructed (Fig. S3a).Likely due to incompatibility with 3-RBD "down" and 1 RBD "up" S-trimer conformations, at least 2 RBDs are observed to adopt "up" position when R1-26 is bound to a S-trimer and only 3:2, 3:3 S-protomer:Fab structures are observed when R1-26 Fab was added in excess (Fig. 2a).The detailed epitope structure also confirmed that the R1-26 epitope does not overlap with the ACE2 binding interface (Fig. S3f, h), indicating that R1-26 Fab is a class 4 antibody (Fig. S3f-h).CR3022 is a well-studied class 4 mAb which cross-reacts to RBDs of SARS-CoV-1 and SARS-CoV-2 39 .We found that although the epitope of R1-26 appears to largely overlap with that of CR3022, R1-26 features a different approach angle towards the RBD compared with CR3022 (Fig. S3a, g, h).Notably, the canonical class 4 antibody CR3022 does not inhibit ACE2 binding and is only weakly neutralizing 16,40 .In contrast, strong ACE2 binding inhibition is observed for R1-26 (Fig. 1b) and other class 4 antibodies such as C118 and C022 40 , as well as S2A4 and S2X35 41 .By modeling, we found that the approach angle of R1-26 is more tilted towards the modeled ACE2 bound to RBD (Fig. S4a), while the approach angle of CR3022 is more tilted away from it (Fig. S4b).In addition, R1-26 and CR3022 adopt different orientations when bound to RBD (Fig. S4a, b).Due to the differences in RBD binding, the simultaneous binding of R1-26 and ACE2 to the RBD would result in steric clashes between R1-26 and the glycan chain attached to ACE2 residue N322 (Fig. S4a).By contrast, CR3022 and ACE2 are able to bind the RBD simultaneously without clashing (Fig. S4b).This is consistent with a previous report indicating that CR3022 does not block ACE2 binding 42 .Differences in RBD binding between R1-26 and CR3022 likely confer superior ACE2 blocking activity to R1-26.A few studies have also linked the neutralization activities of class 4 antibodies (including VL6-57 antibodies S2A4) to their ability to block ACE2 binding 40,41 .
We further found that fusogenic competent native S-R trimer (without proline stabilization) incubated with R1-26 undergoes structural transition into a post-fusion conformation (Fig. 2c).This analysis shows that R1-26 possesses an activity to trigger fusogenic conformational change similar to class 1 mAbs including B38 16 (Fig. 2c) and S230 43 .Interestingly, the class 4 antibody CR3022 does not have this triggering activity (Fig. 2c).Both ACE2 binding inhibition and premature triggering of spike fusogenic change by R1-26 could lead to the inhibition effect observed in a S-protein-ACE2 interaction mediated cell-cell membrane fusion assay by which cell entry inhibition activity The genetic and functional properties of the five isolated VL6-57 light-chain utilizing antibodies, as previously reported 16 .Neutralization titers (IC 50 ) of antibodies were determined against a pseudovirus displaying WT SARS-CoV-2 S-protein.Binding affinities (K D ) of antibodies to SARS-CoV-2 WT RBD were determined by biolayer interferometry (BLI) assays.Germline usage and somatic hypermutation (SHM) analysis of antibodies were performed using IMGT/V-QUEST.b Pairwise binding competition to SARS-CoV-2 RBD was assessed by BLI.A biosensor immobilized with WT SARS-CoV-2 RBD was first saturated with one of the VL6-57 utilizing mAbs until the dashed line before submerging into a solution of another mAb or ACE2 to assess competition.An ACE2 competing antibody B38 (purple) and two non-ACE2 competing antibodies recognizing different RBD epitopes, CR3022 (black) 39 and R1-32 (green) 16 were used as controls.c, d BLI binding curves of R1-26 to 2-fold serially diluted RBD (c) or S-protein (d) solutions of WT SARS-CoV-2 and VOCs.The black lines represent the experimentally recorded sensorgram traces, the colored lines represent corresponding fits.Detailed binding kinetics parameters are summarized in Supplementary Table 1.e Neutralization activities of R1-26 towards SARS-CoV-2 WT and VOC authentic viruses in cell culture (data are presented as mean values ± SD, n = 3).Source data for e are provided as a Source Data file.
of antibody can be assessed 44 (Fig. 2d).Both activities likely contribute to the neutralization activity of R1-26.

VL6-57 light chain can pair with diverse heavy chains to bind a convergent epitope
To further understand antigen binding by VL6-57 light-chain utilizing mAbs, we surveyed the Protein Data Bank (PDB) for structurally characterized VL6-57 light-chain utilizing mAbs bound to SARS-CoV-2 S-protein or RBD.Among a dataset of 376 SARS-CoV-2 S-specific mAb complex structures, which are all experimentally determined by X-ray crystallography or cryo-EM, 12 mAbs utilize VL6-57 encoded light chains, of note, among them, 10 belong to class 4 antibodies (see Fig. S5 for a structure gallery of class 4 antibodies), these account for more than a fourth of class 4 antibodies (n = 35) found amongst the 376 structurally characterized SARS-CoV-2 S-specific mAbs (Fig. S5).Further, germline usage analysis shows that VL6-57 is the most frequently used light-chain germline gene among the class 4 antibodies (Fig. S6a, left panel).An enrichment of VL6-57 is also observed for the F2 group antibodies (equivalent to class 4 antibodies, as defined by deep mutational scanning) 7 (Fig. S6a, right panel).We did not find evidence of preferential VL6-57 gene usage in the global B cell repertoire data of SARS-CoV-2-naive individuals 45 .The above analyses strongly suggest that VL6-57 gene preferentially generates S-specific class 4 mAbs, leading to the hypothesis that VL6-57 light chains may play an important role in antigen binding.
We analyzed buried surface area (BSA) to further understand the contribution of VL6-57 light chains in epitope binding.Among the 376 structurally characterized SARS-CoV-2 S-specific mAbs, heavy  16,43,103 .d R1-26 inhibits spike-ACE2 interaction mediated cell-cell membrane fusion (data are presented as mean values ± SD, n = 3).Two-tailed paired student's t test was performed.Source data for c and d are provided as a Source Data file.
A 3.6 Å resolution structure of the H18 binding interface is derived from a focused refinement of the S1:H18:R1-32 Fab complex (Figs.5f,  S2d).Although the heavy-chain genes are different between H18 (VH4-61) and R1-26 (VH3-7), H18 binds the identified convergent epitope with HCDR3 residues W102 and R104 via hydrophobic contact and cation-π interaction respectively.Consistent with our expectation, the R104 residue within the "WLRG" motif is interacting with Y94 within the LCDR3 "QSYDSS" motif via a cation-π interaction.Antigen binding by H18 HCDR3 is further stabilized by hydrogen bonds from LCDR1 similar to R1-26 (Figs. 5f, S2d).These observations confirm that VL6-57 light chain can pair with multiple heavy chains to target the identified convergent epitope on SARS-CoV-2 spike RBD.

Molecular basis of evasion from the shared VL6-57 antibody response by SARS-CoV-2 Omicron variants
During the structural analysis of VL6-57 antibodies we noticed that residues S371, S373 and S375 are within the convergent epitope in the ancestral SARS-CoV-2 spike RBD (Fig. 3), these residues have been considered as a feature of Omicron variants and shown to alter virological behavior of the virus 28,29 .We also noted that among the 26 RBD substitutions observed in SARS-CoV-2 so far (Fig. 6a), 18 substitution positions are shared among SARS-related-CoVs (Fig. 6b) and among the other 8 substituted positions, S371, S375, T376, and R408 fall within the epitope of the VL6-57 class 4 public antibodies (Fig. 6a, b).Mapping of substituted RBD positions and VL6-57 antibody epitope areas indicates that the epitope is relatively conserved among SARS-related-CoVs.Although R1-26 and H18 could bind and neutralize SARS-CoV-2 WT and early VOCs, both of them fully lose binding and neutralization abilities towards Omicron BA.1 (Figs. 1d, e, 5c).BLI assays show that simultaneous rescue mutations at 371, 373 and 375 on Omicron BA.1 RBD (Omicron-(L371S + P373S + F375S)) are able to completely recover binding by VL6-57 mAbs R1-26 and H18 (Fig. 6c, Supplementary Table 4).This result indicates that the evasion of VL6-57 mAb neutralization is specifically mediated by substitutions at S371, S373, and S375.BLI assays using a series of RBD mutants revealed that single S371L, S373P or S375F mutation greatly reduces binding by R1-26 and H18.Combinations of any two simultaneous mutations (S373P/ S375F, S371L/S375F, S371L/S373P) are able to almost completely abolish R1-26 and H18 binding (Fig. 6c, Supplementary Table 4).Structural comparison of WT, Omicron BA.1 and BA.2 RBDs reveals that backbone conformations at residues 373 to 375 are affected by the S373P introduced proline and the phenylalanine sidechain introduced by S375F (Fig. 6d), these epitope structural changes likely abolish binding by the VL6-57 mAbs.

Discussion
We identify a class of VL6-57 light chain utilizing antibodies with the ability to pair with diverse heavy chains to target a convergent epitope defined by featured Omicron mutations -S371L/F, S373P, S375F.The identified epitope is cryptic in a "down" RBD and strongly hydrophobic, consistent with that HCDR3s of VL6-57 mAbs share hydrophobic residues to interact with the epitope.Notably, within the epitope, residues S371-S373-S375 are located at the entrance to the identified fatty acid binding pocket of sarbecovirus spike RBD.Binding of linoleic acid within the pocket has been associated with "locked" spike conformation with potential functions in virus assembly [36][37][38]57,58 . Beinga highly dynamic protein, SARS-CoV-2 S-protein has been observed to adopt RBD "up" and "down" conformations.D614G substitution became fixed in the spike protein shortly after the SARS-CoV-2 pandemic 36,[59][60][61] , it has been shown that this substitution increased spike stability and a shift towards a more open (RBD "up") S-trimer has been observed by multiple studies 36,59,62 .We speculate that the change in spike dynamics may have increased immune pressure posed by mAbs targeting cryptic epitopes.Indeed, we find that binding of R1-26, H18 and possibly other VL6-57 mAbs to the convergent cryptic epitope not only blocks ACE2 binding 41 , but also prematurely triggers spike conformational change resulting in spike inactivation.When paired with other RBD-specific antibodies, binding of antibodies to cryptic epitopes may also result in spike inactivation by promoting S-trimer disassembly 16,56 .The increased immune pressure may promote the introduction of S371L/F-S373P-S375F in Omicron variants.The S371L/F-S373P-S375F mutations have also been found to stabilize the interaction between up-down RBDs within the Omicron S-trimer, likely further modulating S-trimer dynamics, stability and receptor binding 56,[63][64][65][66] .The identified cryptic epitope of VL6-57 mAbs is relatively conserved among S-proteins of sarbecoviruses 26 .Before the emergence of Omicron variants, the VL6-57 antibody epitope remains  3. c Neutralization activities (data are presented as mean values ± SD, n = 3) of mAb H18 towards SARS-CoV-2 authentic viruses in cell culture.d Structures of S-GSAS/6 P S-trimers in complex with H18 Fabs in different stoichiometries and conformations. e Structres of S-GSAS/6 P S-trimers or S1s in complex with H18 Fabs and R1-32 Fabs.Structures in d, e are low-pass filtered to 12 Å to reveal flexible regions (also see Fig. S2).H18 Fab, NTD, and RBD are highlighted in magenta, blue, cyan, respectively; other structures are colored gray.f Detailed H18 epitope structure. H18-H an H18-L chains are colored in magenta and purple; CDR loops are indicated, selected interacting residues between RBD and H18 are shown and indicated; thick and thin dashed lines indicate cation-π interactions and hydrogen bonds.Source data for c are provided as a Source Data file.
unchanged in early SARS-CoV-2 VOCs or VOIs.We and others also found that some of the VL6-57 mAbs including H18 and many others 67 are able to cross-react with SARS-CoV-1 RBD, confirming that the VL6-57 epitope is relatively conserved among sarbecoviruses and suggesting that cross-reactive VL6-57 mAbs may be used to combat potential emergent sarbecoviruses.On the other hand, the introduction of S371L/F-S373P-S375F demonstrates an extraordinary adaptability of SARS-CoV-2, providing a piece of evidence to suggest that unusual changes of a pandemic virus may occur in response to potential strong population immune pressures including those posed by population antibodies.Studies have shown that single mutation of S371F/L, S373P or S375F strongly impairs spike processing and spike   mediated virus entry 28 .Among them, S375F has the most detrimental effect 29 , and notably, mutation at this position has not been observed among animal sarbecovirus S-proteins 26 (Fig. 6b).
The above observations highlight the likely necessity of the S371-S373-S375 mutations for SARS-CoV-2 Omicron variants, even at the expense of drastic S-protein phenotypic change with associated impact on virus behavior.It has been shown that antigenic changes can alter SARS-CoV-2 receptor binding 68 , for other viruses, antigenic evolution of the pandemic H3N2 virus was found to substantially weaken hemagglutinin (HA) sialic acid binding affecting virological phenotypes 69 .While our study provides evidence to show that the S371-S373-S375 mutations enable SARS-CoV-2 to evade the identified class of VL6-57 public antibodies, other functions for these mutations continue to emerge.A recent study reports that S375F promotes immune evasion by enhancing spike binding to Siglec-9 on macrophages to impair antigen phagocytosis and presentation 30 .A recent preprint also shows that when G339H, S371F, S373P, and S375F mutations are reverted in a reverse-genetic authentic virus carrying Omicron BA.1 S-protein in an ancestral virus genetic background, the Omicron nasal epithelium cell tropism is almost completely abrogated 27 .However, a complex epistatic relationship between S-protein mutations has been suggested-S371L/F, S373P and S375F may only be adaptive when they co-occur with the other Omicronspecific S-protein substitutions 26,27 .Indeed, the correlation between SARS-CoV-2 tissue tropism, pathogenicity, immune evasion, and transmissibility is highly complex and remains ill-defined.Other factors affecting SARS-CoV-2 infection continue to emerge, utilization of membrane-type matrix metalloproteinases (MT-MMPs) and ADAM (a disintegrin and metalloproteinase) has recently been proposed to correlate with nasal cell tropism [70][71][72] .Other Omicron tissue tropism and pathogenesis determinants related to host interferon responses have been recently implicated [72][73][74][75][76][77][78] .Some Omicron attenuation determinants have been mapped outside of the S-protein gene 79,80 .Nevertheless, multiple functions identified so far for the S371L/F, S373P and S375F mutations illustrate vividly that antigenic change of viral surface protein can have complex consequences and a delicate balance among tissue tropism, pathogenicity, immune evasion, and transmissibility is likely pursued by circulating viruses 1 .
It has been reported that SARS-CoV-2 infections induce several classes of germline antibodies of specific VH genes.Due to similarity of equivalent germline antibody genes among humans, induced classes of germline antibodies are widely present in the population.Certain germline antibodies of the same VH gene can bind a convergent spike epitope with similar molecular interactions mediated by antibody residues encoded by germline sequences 45,49,81,82 .Similar phenomena have also been observed for mAbs induced by infections of other viruses [83][84][85] .A common feature of these germline antibodies is that essential paratope interactions are primarily mediated by antibody heavy chains, such feature has been observed for VH3-53 14 , VH1-58 86 , VH1-2 19,20 , VH1-69 16 , and VH1-24 87 SARS-CoV-2 S-specific antibodies.These heavy-chain mediated shared antibody responses have been thought to be a key driving force of SARS-CoV-2 antigenic drift 18,88,89 .In this study, by contrast, the identified population antibody response is mediated by antibodies with shared usage of VL6-57 light-chain gene.Generally, heavy chains dominate interactions in antibody-antigen complexes, while light chains help modulate heavy-chain conformation 90 .It has been shown that antibodies of VH1-69 heavy chains can pair with irrelevant light chains to bind influenza HAs 91,92 .Our structural study reveals that interactions between the VL6-57encoded mAbs and RBD are mainly mediated by germline gene encoded residues in HCDR3 and LCDR1, while germline gene encoded LCDR3 residues support HCDR3 in antigen binding.Although VL6-57 is most frequently observed to pair with VH3-7, VH4-39, and VH4-59 we provide evidence that many other heavy chains can also pair with VL6-57 to bind SARS-CoV-2 RBD.Therefore, VL6-57 germline sequence provides an efficient framework to allow pairing with diverse heavy chains to generate SARS-CoV-2 RBD-targeting antibodies, such an unique property of VL6-57 gene increases the probability of VL6-57 class of public antibodies being induced across a broad population.This study provides a rare example that a light-chain paired with common binding motifs in various heavy chains can largely dictate the binding mode of a class of antiviral antibodies and provides a molecular mechanism to explain the increased identification of VL6-57 Sspecific mAbs among SARS-CoV-2 infected individuals 45,47,48 .
Given the epitope of the VL6-57 mAbs is largely conserved among animal sarbecoviruses, we speculate that the versatility of VL6-57 in generating SARS-CoV-2 RBD antibodies may be unique within the human population.Wide induction of this class of antibodies contribute to selection pressure on SARS-CoV-2 within the human population, but similar response may be absent in other species.Therefore, S371L/F-S373P-S375F represents a specific constellation of immune escape mutations for SARS-CoV-2 adaptation and transmission within the human population.However, it is worth noting that there is still VL6-57 mAb that is resistant to Omicron variants 93 , implying adaptation of antibody immunity to viral mutations.Therefore, exploring the co-evolution of virus and antibody immunity in future works should be beneficial to understanding co-adaptability between virus and host 94 .Overall, study provides a fresh of mediated populational antibody immune pressure against SARS-CoV-2 at molecular level.Our findings lend support to the hypothesis that convergent antibody responses may influence viral antigenic drift.However, establishing a direct causal link between immune pressure and the emergence of specific mutations such as those seen in Omicron variants likely requires additional longitudinal and functional studies.Also, our findings highlight the importance of monitoring viral sequences for mutations in key epitopic regions.Future vaccine strategies could benefit from incorporating such surveillance data to preemptively adjust vaccine compositions in response to emerging escape variants, ensuring continued vaccine effectiveness.

Expression of monoclonal antibody
The antibody heavy-and light-V genes (VH/VL) were cloned into human IgG1 expression vectors using Clone Express II One Step Cloning Kit (Vazyme, China, Cat# C112-02).When the density of HEK293F cells reached 1 × 10 6 cells/mL, equal amounts of heavy-and light-chain plasmids were transfected into HEK293F cells using EZ cell transfection reagent (Life-iLab Biotech, China, Cat# AC04L092).Following transfection, HEK293F cells were cultured in CD 293 TGE medium (ACRO, China, Cat# CM-1156) containing 10% CD Feed X supplement (ACRO, China, Cat# CF-1116-12) at 37 °C in a humidified atmosphere with 5% CO 2 and shaking at 120 rpm.6 days post transfection, supernatants were harvested and clarified by centrifugation.Supernatants were filtered through 0.22-µm filters (Merck Millipore, Germany, Cat# SLGP033NS) before incubation with Protein A Resin (Genscript, China, Cat# L00210) at room temperature for 2 h for antibody affinity purification.After washing, antibodies were eluted from the Protein A Resin using 0.1 M Na-citrate (pH 3.25) and eluents were neutralized immediately with 1 M Tris-HCl (pH 8.8).Antibodies were concentrated in PBS using Amicon Ultrafilter devices (Merck Millipore, USA, Cat# UFC810096) and stored at −80 °C.

Inhibition of SARS-CoV-2 S-mediated cell-cell fusion
As previously described 44 , we utilized a real-time, quantitative live cell split-GFP fluorescence complementation phenotypic assay to analyze the fusion inhibition activity of antibodies.Briefly, effector cells (HEK293T cells expressing SARS-CoV-2 S and GFP1-10 proteins) were seeded into a 96-well-plate (5 × 10 4 cells per well in DMEM).Each well was cultured in the presence or absence of a test mAb at the indicated concentrations for 2 h at 37 °C.Subsequently, target cells (hACE2-293T cells expressing GFP11 protein, 5 × 10 4 cells per well in DMEM) were added uniformly to each well, and fluorescence images were taken 2 h after the addition of target cells using a Nikon fluorescence microscope or a Keyence BZ-X800 all-in-one fluorescence microscope.The GFP area in each well was quantified on Image J, and percentage inhibition of cell-cell fusion was calculated using the following formula: (1 − (E − N)/(P − N)) × 100%."E" represents the GFP area in the experimental group."P" represents the GFP area in the positive control group, where no mAb was added."N" is the GFP area in the negative control group, where effector HEK293T cells only express GFP1-10.Samples were tested in triplicate, and experiments were repeated at least twice.

Protein expression and purification
The extracellular domain of SARS-CoV-2 spike (S) protein (residues 14-1211) with an N-terminal mu-phosphatase signal peptide, a "R" substitution at the multibasic furin cleavage site (R682-R685), and a C-terminal TEV protease cleavage site followed by a T4 fibritin trimerization motif and a His 6 tag was cloned into the mammalian expression vector pCDNA3.1(+),named "S-R", as previously described 43 .To generate the stabilized S-protein for cryo-EM sample preparation, the sequence of S-protein was modified with six prolines at residues 817, 892, 899, 942, 986, and 987 and the furin cleavage site was changed to "GSAS" 95 , named "S-GSAS/6 P".To express the SARS-CoV-2 RBD, residues 319-541 of S-protein were fused with an N-terminal mu-phosphatase signal peptide and a C-terminal 6×His tag before the sequence was inserted into the pCDNA3.1(+)vector.S-protein or RBD expression vector was transiently transfected into Expi293F using polyethylenimine.Proteins were purified using IMAC (immobilized metal affinity chromatography) following previously described protocols 16,37,56 .All proteins were aliquoted, flash-frozen in liquid nitrogen and stored at −80 °C.

Biolayer interferometry
Binding assays were carried out on an Octet RED96 instrument (Sartorius) using a previously established protocol 16 .Briefly, each Protein A biosensors (Sartorius) was pre-equilibrated in PBST buffer (PBS, pH 7.4, 0.02% Tween-20, 1 mg/ml bovine serum albumin) for 10 min.Subsequently, IgG at 11 μg/ml was loaded onto each biosensor to a response level between 1.6-1.8nm.The IgG immobilized biosensors were submerged into twofold serially diluted (200-3.125nM) RBD or S-protein solutions for 300 s to record association.The biosensors were subsequently submerged into PBST buffer for 600 s to record dissociation.For IgG binding to the generated RBD mutants, biosensors immobilized with IgGs were monitored for association (300 s) in RBD solutions at a fixed concentration of 200 nM before the sensors were submerged into PBST to monitor dissociation (600 s).IgGimmobilized sensors were also submerged into a PBST buffer to record references.Data were reference-subtracted and analyzed using Data Analysis HT v12.0.2.59 software (Sartorius) with a 1:1 fitting model for binding to RBDs and 2:1 fitting model for binding to S-trimers.Raw data and fits were plotted in GraphPad Prism v8.0.
Competition assays were performed on a Gator label-free bioanalysis system (GatorBio).2 μg/ml of SARS-CoV-2 RBD (Sino Biological, Cat# 40592-V08B) was immobilized onto the pre-equilibrated Anti-His biosensors (GatorBio).Biosensors were statured with the first antibody for 300 s before submerging into the second antibody or ACE2 solutions for 200 s.Data were analyzed by the Gator data analysis software (GatorBio) and plotted in GraphPad Prism v8.0.
Ligand-induced conformational change of spike protein S-R diluted to 1 mg/ml (7.09 μM) was incubated with ACE2-Fc or antibodies in IgG or Fab form at a 1:1.1 molar ratio at room temperature for 1 h.The samples were subsequently treated with 50 μg/ml proteinase K at 4 °C for 30 min.Each sample was boiled in 5 × non-reducing SDS loading buffer at 98 °C for 5 min to stop the reaction.Samples were separated by SDS-PAGE on 4-12% gradient gels, before protein bands were transferred onto a polyvinylidene difluoride membrane using a semi-dry blotting system.The membrane was blocked with 5% milk in PBST before the membrane was incubated with a primary antibody (rabbit anti-SARS-CoV-2 S2 polyclonal antibody, Sino Biological, Cat# 40590-T62, 1:2500 dilution) in PBST.After extensive washing, the membrane was incubated with the secondary antibody (horseradish peroxidase-conjugated goat anti-rabbit IgG, Beyotime, Cat# A0208, 1:2500 dilution) in PBST.Finally, protein blots were visualized by chemiluminescence using a Pierce ECL Western Blotting Substrate (Thermo Fisher Scientific, Cat# 32106).

Cryo-EM sample preparation and data collection
To generate S-GSAS/6 P:R1-26 Fab, S-GSAS/6 P:H18 Fab complexes, S-GSAS/6 P at 4.43 mg/ml was incubated with R1-26 Fab or H18 Fab at a 1:1 molar ratio.To generate S-GSAS/6 P:H18 Fab:R1-32 Fab complex, S-GSAS/6 P at 4.43 mg/ml was incubated with H18 and R1-32 Fabs at a 1:1:1 molar ratio.After 1 min incubation at room temperature, each 3 μl sample was supplemented with 0.1% octyl-glucoside (Sigma-Aldrich, Cat# V900365) before it was applied onto a 300-mesh holey carboncoated copper grid (Quantifoil, Cu R1.2/R1.3)pre-treated by glowdischarging at 15 mA for 30 s.Each grid was blotted for 2.5 s with a blot force of 4 at 22 °C and 100% humidity before plunge-freezing in liquid ethane using a Vitrobot Mark IV (Thermo Fisher Scientific).The S-GSAS/6 P:R1-26 complex cryo-grid was imaged in a Titan Krios electron microscope (Thermo Fisher Scientific) operating at 300 kV and equipped with Gatan BioQuantum energy filter (slit width 20 eV) and Post-GIF Gatan K3 Summit direct detection camera.Movie stacks were automatically recorded using EPU at a nominal magnification of ×81,000 in super-resolution mode with a calibrated pixel size of 0.5475 Å and nominal defocus values ranged between −0.8 to −2.0 μm.Each stack was fractionated into 38 frames and exposed at a dose rate of 25 e − /pixel/s for 2.4 s resulting in a total dose of ~50 e − /Å 2 .The S-GSAS/6 P:H18 complex cryo-grid was imaged in a Talos Arctica electron microscope (Thermo Fisher Scientific) operating at 200 kV.Using the SerialEM v3.8.7 software, movie stacks were recorded at a nominal magnification of ×45,000 on a K3 direct detection camera (Gatan) in super-resolution mode with a calibrated pixel size of 0.44 Å with nominal defocus values ranged between −0.8 to −2.5 μm.Each stack was fractionated into 27 frames and exposed at a dose rate of 24.4e − / pixel/s for 1.89 s, resulting in a total dose of ~60 e − /Å 2 .The S-GSAS/ 6 P:H18:R1-32 complex cryo-grid was imaged in a 300 kV Titan Krios electron microscope (Thermo Fisher Scientific) equipped with a SelectrisX energy filter (slit width 10 eV) and a Falcon 4 direct electron detector.Movie stacks were automated collected using EPU software with the electron event representation (EER) mode at a nominal magnification of ×130,000 with a pixel size of 0.93 Å and nominal defocus values ranged between −0.6 to −2.0 μm.Each stack was recorded and exposed at a dose rate of 7.51 e − /pixel/s for 5.79 s resulting a total dose of ~50 e − /Å 2 .These settings yielded EER stacks consisting of 199 frames.All movie stacks were imported into cryoS-PARC live (v3.3.2) 96 for pre-processing, which includes patched motion correction, contrast transfer function (CTF) estimation and bad images rejection.Movie stacks of S-GSAS/6 P:R1-26 complex and S-GSAS/ 6 P:H18 complex datasets were binned 2× resulting in pixel sizes of 1.095 Å and 0.88 Å, respectively.

Cryo-EM data processing
Data processing was carried out using cryoSPARC v3.3.2/v4.2.0.After bad image removal, particles were picked by blob-picking on 4343 S-GSAS/6 P:R1-26, 4822 S-GSAS/6 P:H18 and 9861 S-GSAS/6 P:H18:R1-32 images.After 2D Classification, several good 2D classes were selected as the templates for template-picking, resulting in initial datasets of 2,359,986 particles for S-GSAS/6 P:R1-26; 2,357,740 particles for S-GSAS/6 P:H18; and 3,335,676 particles for S-GSAS/6 P:H18:R1-32.Template-picked particles were extracted for two or more rounds of 2D Classification to remove contaminants and low-quality particles.For S-GSAS/6 P:H18:R1-32, particles with intact S-trimers were reextracted with a larger box size before the last round of 2D Classification.After several rounds of 2D Classification, well-aligned particles with intact S-protein features were selected and subjected to Ab-initio Reconstruction to generate 3D models.Particles of initial models with complete S2 but incomplete S1 were further classified using Ab-initio Reconstruction with class similarity value 0.3 for the S-GSAS/6 P:R1-26 dataset, 0.7 for the S-GSAS/6 P:H18 and S-GSAS/6 P:H18:R1-32 datasets to separate different conformations.
3:2 and 3:3 S-protomer:Fab structures were observed in both S-GSAS/6 P:R1-26 and S-GSAS/6 P:H18 datasets.Head-to-head aggregates of 2 S-trimers bound by 3 Fabs, giving S-protomer:Fab stoichiometry of 6:6, were observed in 2D Classifications and Ab-Initio Reconstructions of both S-GSAS/6 P:R1-26 and S-GSAS/6 P:H18 datasets.Particles of aggregate classes were selected and combined before particles were re-extracted with a larger box size.A S-trimer aggregate structure from the S-GSAS/6 P:R1-26 dataset was obtained by a round of 2D Classification followed by Ab-Initio Reconstruction.Similarly, 2 S-trimer aggregate structures were obtained from the S-GSAS/6 P:H18 dataset.To obtain a higher resolution structure of the interface between RBD and R1-26 Fab, particles of the 3:3 S-GSAS/6 P:R1-26 structure were refined again with applied C3 symmetry by Non-uniform Refinement.Refined particles were subjected to Symmetry Expansion and Density Subtraction before a Local Refinement was carried out focusing on a region containing RBD and Fab-V.
All maps with resolutions higher than 4 Å were processed by a second round of Non-uniform Refinement with defocus refinement and global CTF refinement to improve map quality.All resolutions were estimated at the 0.143 criterion of the phase-randomizationcorrected Fourier shell correlation (FSC) curve calculated between two independently refined half-maps multiplied by a soft-edged solvent mask in RELION v4.0 96 .Additional data processing details are summarized in Figs.S1 and S2, and Supplementary Table 2.

Model building and analysis
A previously determined structure of SARS-CoV-2 S-trimer in complex with 3 R1-32 Fabs and 3 ACE2s (PDB: 7YEG) 16 or a structure of SARS-CoV-2 S1 in complex with a YB9-258 Fab and an R1-32 Fab (PDB: 8HC5) 56 was used as starting model.An antibody light chain (PDB: 7D6I) was used as the starting model for R1-26 and H18 Fab light chains.Starting models for R1-26 and H18 Fab heavy chains were generated from heavy-chain structures in PDBs -5X8M and 7VSU respectively.Starting models were fitted into final refined maps in UCSF Chimera v1.14 97 .Iterative model building and real space refinement were carried out in Coot v0.9.6 98 and PHENIX v1.20.1-4487 99.Model refinement statistics are summarized in Supplementary Table 2. Interface analyses were performed in QtPISA v2.1.0 100.Structure figures were generated in UCSF Chimera v1.14.

Analysis of the structurally characterized SARS-CoV-2 S-specific mAbs
A total of 376 antibody:RBD/S-protein complex structures were curated from the PDB (https://www.rcsb.org/)before May 18, 2023.LCDR2 loop structures of 2-7 (PDB 7T3M, EMD-25663), 553-15 (PDB 7WO7, EMD-32641) were lightly adjusted in Coot v0.9.6 98 according to the corresponding density maps and refined in PHENIX v1.20.1-4487 99.Epitope residues, paratope residues, and buried surface area (BSA) were determined or calculated using the PDBe PISA server (https:// www.ebi.ac.uk/msd-srv/prot_int/).BSA for each epitope residue is assigned as a feature of a certain antibody and used to construct a feature matrix M AxB for downstream analysis, where A is the number of antibodies and B is the number of features (amino-acid length of spike: 1273).Therefore, a 376 × 1273 BSA matrix was obtained, which was subsequently used as input for epitope classification with the R package UMAP (v0.2.9.0).Uniform Manifold Approximation and Projection (UMAP) is an algorithm for dimensional reduction.After the reduction analysis, clustering was performed using the K-means algorithm.UMAP and K-means clusterings were conducted independently.Antibodies clustered into the same cluster as the well-documented mAb CR3022 are defined as class 4 antibodies.After this workflow, 35 class 4 mAbs were identified among the 376 structurally characterized S-specific mAbs (Fig. S5).In addition, we found that class 4 epitope is largely overlapping with the F2 epitope which was previously defined by deep mutational scanning 7 .Thus, a total of 69 F2 epitope mAbs were also included for germline gene usage analysis.Germline gene usages of these structurally characterized S-specific mAbs were inferred using IMGT/V-QUEST (https://www.imgt.org/).Visualization of the germline heavy and light gene usage and pairing among the 35 class 4 and the 69 F2 antibodies were performed in R platform v4.2 using the R package circlize v0.4.10 101 .

Bioinformatic analysis
To investigate features of VL6-57 utilizing mAbs, a total of 290 VL6-57 utilizing mAbs were extracted from the CoV-AbDab database 46 for downstream analysis.Analysis of germline gene usage and CDR3 length were perform using the built-in function in R platform v4.2.CDR3 amino-acid composition analysis was performed using the R package ggseqlogo 102 .VL6-57 utilizing mAbs with 12-AA length HCDR3 were used for lineage structure reconstruction.R package igraph was employed for visualization of the lineage structure of VL6-57 utilizing mAbs with 12-AA length HCDR3s.After the above workflow, the results showed that WLRG motif within 12-AA long HCDR3 and QSYDSS motif within LCDR3 are the convergent signatures of the VL6-57 mAbs.To determine the occurrence of VL6-57 mAbs with paired HCDR3 "WLRG" motif, we analyzed nearly 2 billion IgH and IgL sequences from 3 previously described datasets 13,49,50 .These datasets have been deposited by us in the National Genomics Data Center (https://bigd.big.ac.cn/),China National Center for Bioinformation (CNCB) under accession numbers PRJCA003775, PRJCA007067, and PRJCA017560.HCDR3 "WLRG" motif-containing IgH sequences were defined as those sequences which have 12-AA long HCDR3s with a "WLRG" motif.VL6-57 IgL sequences were defined as those sequences that use the germline gene VL6-57 with 9 to 10-AA long LCDR3s, and encoding a "QSYDSS" motif within LCDR3 region.The abundance of "WLRG" motifcontaining IgH sequences or VL6-57 IgL sequences were normalized by counts per million.The clonal expansion of B cells expressing "WLRG" motif-containing IgH or VL6-57 IgL sequences in COVID-19 patients after SARS-CoV-2 infection were tracked.The divergence from germline genes of all IgH or IgL sequences was equal to somatic hypermutations, and sequence identity to queried sequence was calculated using the R package Biostrings v2.60.2 (http://bioconductor. org/packages/release/bioc/html/Biostrings.html).All IgH or IgL sequences were plotted as a function of sequence somatic hypermutations (x axis) and sequence identity (y axis) to heavy or light chains of R1-26, R1-30, R2-3, R2-6, or R2-7 with a color gradient indicating sequence density.To determine the frequency of B cells expressing VL6-57 mAbs with paired HCDR3 "WLRG" motif, we analysis the single-B V(D)J sequences from SARS-CoV-2-exposed and -naive individuals [51][52][53][54][55] and FACS-sorted RBD + B cells from COVID-19 convalescents 56 .These datasets are available from the Gene Expression Omnibus database under accession numbers: GSE230227, GSE171703, GSE158038, GSE158055; http://www.microbiome-bigdata. com/project/SARS-CoV-2/; and National Genomics Data Center (https://bigd.big.ac.cn/) under the accession number: PRJCA012020.B cells expressing VL6-57 mAbs with paired HCDR3 "WLRG" motif were defined as those B cells expressing heavy chains with a 12-AA long HCDR3 containing a "WLRG" motif and VL6-57 light chains with a "QSYDSS" motif in LCDR3.A schematic diagram summarizing the bioinformatic analysis workflow is shown in Fig. 4a, which is created with BioRender.com,released under a Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International license.

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Fig. 1 |
Fig.1| Characteristics of the 5 VL6-57 light-chain utilizing antibodies.a The genetic and functional properties of the five isolated VL6-57 light-chain utilizing antibodies, as previously reported16 .Neutralization titers (IC 50 ) of antibodies were determined against a pseudovirus displaying WT SARS-CoV-2 S-protein.Binding affinities (K D ) of antibodies to SARS-CoV-2 WT RBD were determined by biolayer interferometry (BLI) assays.Germline usage and somatic hypermutation (SHM) analysis of antibodies were performed using IMGT/V-QUEST.b Pairwise binding competition to SARS-CoV-2 RBD was assessed by BLI.A biosensor immobilized with WT SARS-CoV-2 RBD was first saturated with one of the VL6-57 utilizing mAbs until the dashed line before submerging into a solution of another mAb or ACE2 to assess

Fig. 2 |
Fig. 2 | Structural analysis of R1-26 epitope binding and biochemical activities of R1-26.a Structures (low-pass filtered to 12 Å) of S-GSAS/6P S-trimers in complex with R1-26 Fabs in different stoichiometries.R1-26 Fab, NTD and RBD are highlighted in magenta, blue and cyan, respectively; the rest of the S-trimer is colored gray.b Epitope of R1-26 on spike RBD.R1-26-H and R1-26-L chains are colored in magenta and purple.Detailed interactions between R1-26 and RBD are shown in the dashed boxes.CDR loops are indicated, selected interacting residues in antibody-RBD interface are shown.Thick and thin dashed lines indicate cation-π and hydrogen bond interactions; The solid line represents a salt bridge.c Ligandinduced conformational change assays to probe the induction of post-fusion structures, presence of proteinase K resistance core is indicative of post-fusion structure16,43,103 .d R1-26 inhibits spike-ACE2 interaction mediated cell-cell membrane fusion (data are presented as mean values ± SD, n = 3).Two-tailed paired student's t test was performed.Source data for c and d are provided as a Source Data file.

Fig. 3 |
Fig. 3 | Epitopes and binding modes of VL6-57 antibodies.a-i Structures of R1-26 (a), FP-12A (b), IY-2A (c), 553-15 (d), 2-7 (e), S2A4 (f), P5S-3B11 (g), GH12 (h), and 3D11 (i) bound to their epitopes.Fab-VH, Fab-VL, and RBD are colored in magenta, purple, and cyan, respectively.Selected interactions engaged by HCDR3, LCDR1, and LCDR3 are shown.All shown antibodies bind almost the same surface area on RBD, hydrophobic, particularly aromatic amino-acid residues are found in all the HCDR3 loops of the shown VL6-57 antibodies for epitope recognition.All LCDR1 loops have the same interactions with RBD.j Comparison of RBD binding by the shown VL6-57 mAbs and ACE2.k Superposition of HCDR3 and LCDR1 loops of the shown VL6-57 mAbs.Variability is observed for the HCDR3 loops while the LCDR1 loops show high similarity.

Fig. 4 |
Fig. 4 | Analysis of VL6-57 light-chain sequences and associated heavy-chain sequences in antibody repertoires.a A schematic diagram showing the search workflow for VL6-57 light chain and the associated "WLRG" motif containing heavychain sequences from CoV-AbDab or published antibody repertoire datasets.Figure 4a, created with BioRender.com,released under a Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International license.b Heavy-chain sequence analysis of the 290 VL6-57 utilizing mAbs identified from the CoV-AbDab.A histogram showing the HCDR3 length distribution among the 290 mAbs (bottom).Sequence logo plot showing consensus HCDR3 amino acids of the VL6-57 utilizing mAbs with 12-AA long HCDR3s (top left).A histogram showing the heavychain germline gene usage among the VL6-57 utilizing mAbs with 12-AA long HCDR3s (top right).Lineage structure of the VL6-57 utilizing mAbs with 12-AA long HCDR3s based on similarity of HCDR3 sequences (boxed).When the similarity of any two mAbs' HCDR3s is ≥80%, they are linked by a line.mAbs containing the HCDR3 "WLRG" motif are colored red.c Pie chart showing the occurrence of "WLRG" motif containing IgH sequences in a cohort of 33 COVID-19 convalescents (sample size, n = 63) and 24 naive donors (ND, n = 24).d Read-count comparison of "WLRG" motif containing IgH sequences in the IgH repertoires of COVID-19 convalescents (n = 63) and naive donors (n = 24).e Read-count dynamics of "WLRG" motif containing IgH sequences in the IgH repertoires of COVID-19 convalescents, each line connects data of different timepoints from the same patient.f Read-count comparison of "WLRG" motif containing IgH sequences in the IgH repertoires of COVID-19 convalescents and naive donors across different timepoints (n = 24, 32, 31 for ND, 1-13, and 14-28 group, respectively).g Pie chart showing the isotype distribution of "WLRG" motif containing IgH sequences in COVID-19 patients and naive donors.h Read-count comparison of "WLRG" motif containing IgH sequences between COVID-19 convalescents (n = 63) and naive donors (n = 24) across different isotypes.i Read-count comparison of "QSYDSS" containing VL6-57 light-chain sequences in the IgL repertoires of COVID-19 convalescents and naive donors (n = 22, 23 for ND and 1-13 group, respectively).Antibody sequence analyses in c-i used published bulk sequencing datasets13,49,50 .j The occurrence of VL6-57 mAbs with paired heavy chains containing the HCDR3 "WLRG" motif in the single-B V(D)J repertoires of COVID-19 patients, vaccinees and naive donors[51][52][53][54][55][56] .Data are presented as mean values ± SEM in d, f, h, and i.The two-sided chi-square test was performed in c and g.Two-tailed unpaired student's t test was performed in d, f, h, and i. Source data for d-f, h, and i are provided as a Source Data file.

Fig. 5 |
Fig. 5 | Characteristics of recombinant mAbs derived from SARS-CoV-2 naive donors.a WT SARS-CoV-2 RBD binding by five recombinant mAbs were assessed by BLI.Black and colored lines represent experimentally recorded sensorgram traces and corresponding fits.Twofold serially diluted RBD solutions (200 nM to 3.125 nM) were used in the BLI assays.b Binding of mAb H18 to SARS-CoV-1, Pangolin CoV GD1, and Bat CoV RaTG13 RBDs were measured by BLI.Kinetic parameters for a and b are summarized in Supplementary Table3.c Neutralization activities (data are presented as mean values ± SD, n = 3) of mAb H18 towards SARS-CoV-2 authentic viruses in cell culture.d Structures of S-GSAS/6 P S-trimers in

Fig. 6 |
Fig. 6 | variants VL6-57 antibodies by mutations at residues S371-S373-S375.a Substituted residues among SARS-CoV-2 variants are highlighted in red on the surface of SARS-CoV-2 RBD.Among the 26 substitutions, 18 are shared among SARS-related-CoVs and are labeled in red, the other 8 substitutions are labeled in orange.b Substituted residues among RBDs of SARSrelated-CoVs are highlighted in red.Among the 72 substitutions, 18 are shared with SARS-CoV-2 variants and are labeled in red.Epitope of R1-26 on RBD is encircled by cyan outlines in a and b. c Binding of R1-26, H18, and S2A4 to a series of Omicron BA.1 RBD mutants carrying 1-3 rescue mutations at S371, S373, and S375.BLI binding kinetic parameters are summarized in Supplementary Table4.d Comparison of the VL6-57 mAb epitope structures on WT and Omicron BA.1 (left) or BA.2 (right) RBDs.Superposition shows differences in amino-acid sidechain orientations and main chain backbone structures.Source data for a and b are provided as a Source Data file.